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On forced oscillations in groups 
of interacting nonlinear systems 


Ivan Polekhin 


Abstract 

Consider a periodically forced nonlinear system which can be presented as a collection of 
smaller snbsystems with pairwise interactions between them. Each snbsystem is assnmed 
to be a massive point moving with friction on a compact snrface, possibly with a bonndary, 
in an external periodic held. We present snfhcient conditions for the existence of a periodic 
solntion for the whole system. The resnlt is illustrated by a series of examples including a 
chain of strongly coupled pendulums in a periodic held. 

Keywords: periodic solution, Euler-Poincare characteristic, nonlinear system, coupled 
pendulums, nonlinear lattice 


1. Brief introduction 

The phenomenon of forced oscillations has been studied since at least 1922, when G. Hamel 
proved [Ij that the equation describing the motion of a periodically forced pendulum have at 
least one periodic solution. Various results, which generalize and develop P, have appeared 
in the literature since this paper; see, for example, p El SI El E]. Forced oscillations in a 
system of coupled planar pendulums and its generalizations has been studied in PE] under 
the assumption of certain symmetry properties for the forcing terms. The existence and 
multiplicity of periodic solutions for coupled systems also discussed in [9]. 

At the same time, a nonlinear lattice is a cornerstone model in nonlinear physics and 
it is widely used for analytical and computational purposes. See, for example, [iniiiiiEi 
[13]. In our work we present a result which lies in the intersection of the both mentioned 
research areas and can be useful for studying forced oscillations in a nonlinear lattice and 
its generalizations. The result given in the paper continues a previously reported result [TT| . 

In the paper the following systems are considered. Let us have several compact smooth 
manifolds, possibly with boundaries, with non-zero Euler-Poincare characteristics. Suppose 
that for each manifold there is a massive point moving on it with viscous friction. All 
points are in an external periodic held and may interact with each other. It is allowed that 
the interactions may be of different types and also may be arbitrarily strong. We present 
sufficient conditions for the existence of a periodic solution for such a system. 
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We would like to note that when the considered manifolds are closed, it isn’t hard to prove 
that there exists a periodic solution in the system: it directly follows from the Lefschetz-Hopf 
theorem. Yet in applications it appears to be useful to consider manifolds with boundaries 
to prove the existence of a periodic solution and estimate it. One of the examples of such 
kind, which we consider further below, is a widely used in nonlinear studies model of coupled 
oscillators on a line. Therefore, our result generalizes the approach based on an application 
of the Lefschetz-Hopf theorem to the important case of compact manifolds with boundaries. 

Our result is based on a theorem by R. Srzednicki, K. Wojcik, and P. Zgliczyhski [15] 
and provides an illustrative geometrical approach to study periodical oscillations. Since the 
result is presented in a coordinate-free form, it also avoids possible shortcomings of purely 
analytical approaches and can be applied to the systems with complex topology of the phase 
space. 

The main theorem is illustrated by a series of examples from mechanics including a 
system of an arbitrary number of strongly coupled pendulums in a periodic external held. 


2. Main result 

2.1. Governing equations 


The equations introduced in this subsection — which we are going to use further below 
in the paper — generalize the governing equations for a mechanical system consisting of 
massive points moving with friction-like interaction on compact surfaces. For the sake of 
simplicity, we assume that all manifolds and considered functions are smooth (i.e. C°°). 

Let Mi be a compact connected one- or two-dimensional manifold (possibly with a bound¬ 
ary), where i = 1,..., n. For alH = 1,..., n, there is a point moving in an external force held 
on the manifold Mj. We also assume that there is a friction-like force acting on the point. 

In our further consideration, we will study the behavior of the system in vicinities of dMi 
and it will be convenient to consider enlarged manifolds . Let be a boundaryless 
connected manifold such that Mj C , dimMj = dimMj^. Also suppose that every 
manifold is equipped with a Riemannian metric (•, ■)i. 

If the points do not interact, one can consider the following independent equations of 
motion 



fi{t, Qi, Qi) + Qi, (ji), i 


l,...,u. 


( 1 ) 


where V* means the covariant diherentiation with respect to the corresponding metric (■, ■)i 
for the f-th point moving on fi \ M/TZ x TMf — )■ TMf and fi{t,qi,qi) G Tg.M^ for 
any t, qi, qi] M/TZ x TMi^ —)■ TM+ corresponds to a friction-like force. 

As we said before, ([^ includes the case of a mechanical system of massive points in 
an external held. Indeed, in this case (•,-)i is the corresponding kinetic metric and, for 
given t, qi and qi, fi{t,qi,qi) and are the dual vectors to the corresponding 

generalized forces. Note that from ([^ it follows that if there is no forces acting on the Ath 
point then (j^ becomes = 0, which is the equation of the geodesic motion. 
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In general form, when the interactions are took into acconnt, we consider the following 
equations of motion 


= fiit, Qi, Qi) + qi) + /, 


cinteraction 


{t,qi,qi,...,qn,qn), i = (2) 


where X TM^ X ... X TM+ —)■ TM^ and gi, gi,..., gn, <?n) G 

Tq^M^ for any t, g^, g^, j = 1,..., n. Here the functions ^ being dehned on M/TZ x 

TMi X • ■ • X TM^, may describe more general forces than the pairwise interactions usually 
considered in applications. 

As usual, we say that the system has a T-periodic solution if there are functions g^ : M —>■ 
which satisfy @ and qi{t + T) = gi(f), for all t e [0,T), i = 1, ...,n. 

Finally, rewrite (|^ as follows 


Qi = Pi, 


^nPi = Mt, qi,Pi) + gi,pi) + /* 


cinteraction 


(f,gi,Pi, ...,qn,Pn), 


i = (3) 


In the rest of the paper, we will assume that ^interaction gg^i^jgfy following 

conditions: 


(HI) There exist constants d* > 0 such that 

(f frictions. ^ \ 

\ji qi, Pi), Pi/i 


sup 

te[0,T],gieMi 
{Pi ,Pi)i>di 


{Pi,Pi)i 


< 0, i = 1,..., n. 


(4) 


(H2) The functions fi and bounded for all i = 1, ...,n. 

Remark 2.1. If we consider a mechanical system of massive points moving on surfaces in 
with viscous friction with arbitrarily small friction coefficients then (|^ is satisfied. 

Lemma 2.2. Suppose that conditions (HI) and (H2) are satisfied, then for some Ci > 0, 
z = 1, ...,n along the solutions of 


_T- 

dt * 


Ti=Ci 


< 0 , Ti = {pi,pi)i = \\pi 


Proof. By direct calculation from (|^, we have 


d 


jTi = 2{fi{t, qi,pi) + q.^p.) + fi 


^ 2{pi,pi)i 


\\fiit,qi,Pi)\\i ^ \\fi 


interaction 


friction/, \ \ 

{t,qi,Pi),Pi)i 


\\Pi\\i 


it,qi,Pl,- ,qn,Pn)\\i ^ {fi 


\\Pi\\i 


{pi,Pi)i 


Since fi and bounded, then from (|^ we obtain dfi/dt < 0 provided q is large 

enough. □ 


M+ 
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2.2. Auxiliary constructions and results 

The approach developed in [15] is based on the ideas of the Wazewski method [T6| and 
the Lefschetz-Hopf theorem. In this subsection we introduce some dehnitions and a result 
from [TS] which we slightly modify for our use. 

Let uiMxM—j-TMbea time-dependent vector field on a manifold M 

X = v(t, x). (5) 

For to eM. and Xq G M, the map t h-)■ x(t, to, Xq) is the solution for the initial value problem 
for the system ([^, such that x(0, to, xq) = xq. If IF C M x M, f G M, then we denote 

Wt = {xE M: {t,x) G IF}. 

Definition 2.3. Let IF C M x M. Dehne the exit set 1F“ as follows. A point (^ 052 : 0 ) is in 
1F“ if there exists 5 > 0 such that {t + to, x{t, to, xq)) ^ IF for all t G (0, 6). 

Definition 2.4. We call IF C M x M a Wazewski block for the system ([ 5 ]) if IF and 1F“ 
are compact. 

Definition 2.5. A set IF C [a,b] x M is called a simple periodic segment over [a,b] if it 
is a Wazewski block with respect to the system (|^, IF = [a,b] x Z, where Z C M, and 
for any fi, ^2 G [a, b). 

Definition 2.6. Let IF be a simple periodic segment over [a, b]. The set IF = [a, b] x IF^ 
is called the essential exit set for IF. 

In our case, the result from [T5] can be presented as follows. 

Theorem 2.7. JT^ Let W be a simple periodic segment over [a,b]. Then the set 

U = {xo G Wa : x{t — a, a, Xq) G IFi \ W^ for all t G [a, 6]} 

is open in IVa and the set of fixed points of the restriction x{b — a, a, ■)\u'. U —)■ Wa is 

compact. Moreover, if Wa and Wf are ANRs then the fixed point index of x{b — a, a, •)|{/ 

can be calculated by means of the Euler-Poincare characteristic of W and Wf as follows 

md{x{b - a, a, ■)\u) = x(.Wa) - x(lF“). 

In particular, if x(lFa) — x(1FT) 7 ^ 0 then x{b — a, a, ■)\u has a fixed point in Wa. 
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2.3. Main theorem 

In this subsection, we prove our main result and illustrate it with a series of examples. 
Below we use the following notations: M = Mi x ... x and x ... x M+. 

Theorem 2.8. Suppose that for ^ the following conditions are satisfied 

1. The Euler-Poincare characteristic of Mi is non-zero for all i = 1, ...,n. 

2. (HI) and (H2) are satisfied. 

3. For any to G M, {(li,Pi, ■■■, QnHn) ^ T{dM) there is an e > such that 

q{t,to,q(,p[,...,q^,p^J ^ M, for all te{0,e). (6) 


Then there exists a solution {qi,pi, ...,qn,Pn) ■ 1^ —t TM+ of such that for all i = 1, ■■■, n 
Qi{i) = Qi{t+ T), pi{t) = pi{t+ T), qi{t) e Mi\dMi, for all teR. 


Proof. Consider the following compact subset W of [0,T] x TM^ 

W = {0 ^t ^T, {qi,pi, ...,qn,Pn) G TM+: (gi, ...,qn) G MiX...xM„, {pi,Pi)i ^Ci,i = 1, ...,n} 


where q > 0 are the constants obtained from lemma |2.2| 

From lemma 2.2 we also have that if (t,qi,pi, ...,qn,Pn) G W then (gi,...,g„) G dM. 
Since 

dM = dMi X M 2 X ... X MnU ... U Mi x ...M„_i x (9M„, (7) 


then for {t,qi,pi, ...,qn,Pn) G W there exists i such that qt G dMi. Let G Tq.Mf' be 
a normal vector to dMi at point g* G dMi such that, for all pi G Tg^Mf', {i'qi,Pi)i > 0, the 
solution starting from (t, gi,pi, ..., qi,Pi, ..., qn,Pn) at least locally leaves M. From the above 
dehnition of Ug ., we obtain for any j 


W D {0 < f ^ T, (gi,pi, ...,g„,p„) G TM+: 

(qi, ■■■,qn) G M, {pi,Pi)i ^ Ci,i = 1, ...,n,qj G dMj, {i'q^,Pj)j > 0}. 

Let us denote by Id- , j = 1,..., n the following set 

V~~ = {0 ^ t ^ T, (gi,pi, ...,qn,Pn) G TM+: 

(gi, ■■■,qn) G M, {pi,Pi)i ^Ci,i = 1, ...,n, {i'q^,Pj)j > 0, g^ G dMj). 

Since we assume §. then W is compact: 


W— = I/f-U...UI/”. 


( 8 ) 


Therefore, W is also compact and IT is a simple periodic segment over [0,T]. Clearly, 
T- is homotopic to Mi x ... x dMi x ... x M„. Moreover, V, fl... fl K is homotopic to 

J J J\ Jm 

Ml X ... X dMj.^ X ... X dMj^ x ... x M„. 


5 





Figure 1 : For a one-dimensional manifold parametrized by g S [<71,92], bF has the above form. 

If Mi has a nonempty boundary, then dMi consists either of a hnite number of curves 
that are homeomorphic to circles (x(cIMj) = 0) or it is a two-pointed set (x(cIMj) = 2). 
Suppose, without loss of generality, that for some k x{dMi) = 0 iS i > k, i.e. there are 
k ^ 0 one-dimensional manifolds in {Mj}. From ([^ we have 

x(w-)^ Y.x(v--)- Y. x(v--r<v--) + ... + {-ir-\{Vf-n...nv--) 

= E - E n... n y-) 

= x(Mi X ... X Mn) ■ (201 - 2^Cl + ... + 2\-lf-^Cl) = x{M) ■ (1 + (-1)"+'). 
Here we use that if ji < ... < jm ^ k^ then 

x(Mi X ... X dMj^ X ... X dMj^ x ... x M„) = 2’”x(Mi x ... x M^). 

Finally, x(IFo) = x(M) = x{Mi)...x{Mn) + 0, x(IF”) = x(IFo“). Therefore, 

x(W^o)-x(W^o“)^0 


and we can apply theorem 2.7 


Let us now illustrate the theorem on examples. 


□ 


Example 2.9. Consider a hnite number of planar pendulums moving with viscous friction 
in a gravitational held (Fig. [^. Let be the radius-vector from a hxed origin to the massive 
point of the i-th pendulum. Suppose that the pivot points of the pendulums are moving 
along a horizontal line in accordance with a T-periodic law of motion h: M/TZ —)• M, which 
is the same for all pendulums. Let us also assume the following: for any two pendulums 
the distance between their pivot points is greater than the sum of their lengths, i.e. the 
pendulums are disjoint; for any two pendulums there is a repelling force Fjj acting on the 
massive point of the Fth pendulum from the j-th pendulum [Fij is parallel to r* — rj). 
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Figure 2 : When the i-th pendulum is horizontal, the repelling forces acting on it are directed downward. 


The equations of motion can be presented in the following form: 


= Pi, 

h Q 1 V , . 

Pi = -- cos + - sm Lpi - 'yiPi H-- > {Fij, e^J. 

mik ^ 


(9) 


U 


k 




Here m, is the mass of the t-th massive point; U is the length of the Tth pendulum; 7i > 0 
is the viscous friction coefficient for the Tth pendulum (divided by the mass ruj); e^. is the 
following unit vector e<^. = (cosp,, — sinp*); for given ipi, ipj, pi, pj, we assume that 


Fij i}Pi,Pi, Pj, Pj ) 


fij {Pi, Pj), 


where fij ^ 0 is a smooth real-valued function. Then there exists a T-periodic solution of (|^. 

Indeed, let Mi be a compact manifold defined by the inequalities —7r/2 ^ pi ^ 7r/2, 
xiMi) = 1. From (|^ it follows that 

1 


Pi\,„.-^/o „.-n 1 + j ®V’i)|^.=vr/2 ^ 


Ii5i=7r/2,p,=0 1 . 


1 


PA<pi=-w/2,pi=0 
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^ m./. 


( 10 ) 


(/Ji = - 7 r /2 


< 0 . 




Therefore, theorem 2.8 can be applied. Let us also note that along the obtained periodic 
solution we always have pi G (—7r/2,7r/2), i.e. the pendulums always remain above the 
horizontal line. 


Remark 2.10. Note that if in the previous example we replace the half-circles by any compact 
manifolds with boundaries such that their Euler-Poincare characteristics are non-zero, the 
boundaries are in a horizontal plane, the rest of the manifolds are above the horizontal plane 
and they are vertical at the boundaries, then, for such a system, there also exists a periodic 
solution. Moreover, we can consider arbitrary external periodic helds acting on the massive 
points instead of the periodic motion of the surfaces. If the external helds are directed 
downward at the horizontal plane then there exists a periodic solution (Fig. [^. 
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Figure 3 : The points are repelling; the external fields Fi and F2 are periodic and directed downward at the 
boundaries. 


Remark 2.11. One can also consider systems of conpled pendnlnms in which the pendnlnms 
attract each other. An example of snch a system is presented in the fignre below (Fig. |^. 
Here we assume that the massive points of the pendulums are charged, they interact with 
the Coulomb interaction and the pendulums are located in a parallel electric field E. This 
system is similar to example 2.9 and also has a periodic solution. 



Figure 4 : An example of the system with a periodic solution. The pendulums attracts each other and their 
pivot points are moving periodically along the horizontal line. 


Example 2.12. Consider a chain of coupled oscillators in an external periodic field. We 
suppose that the oscillators are located along a straight line, i.e. every oscillator has at 
most two neighbours. For the boundary oscillators, which have only one neighbour, we 
assume that the corresponding massive points are fixed. We also assume that the oscillators 
are moving with viscous friction and the interaction between them is described by a Morse 
potential (Fig.[^. For simplicity, let = 1 and 7 > 0 is the same viscous friction coefficient 
for all oscillators. The equations of motion can be written as follows: 


dV{xi-Xi_i) dV{xi+i-Xi) , , 

Xi = -jXi --- \-F{t,Xi), z = l,...,n, 


dxi 


dxi 


( 11 ) 


where 


V{x) = 
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F: M/TZ X M —)■ M and 7 > 0, 5 > 0, Xq = 0 and Xn+i = (2n + 1)((5 + a), a ^ ln2 is a 
parameter. 


Ml 


M 2 

X3 

M3 







6 + a 

5 CL 


Figure 5: For n = 3 we have three compact manifolds with boundaries. The manifolds Mi are distanced 
from each other so the oscillators attract one another; Xi is directed outward at the boundary dMi. 


Let us now show that if the function F satishes the following condition 

(_l)^+i . for all t G M/TZ, k = 1,2n, 


( 12 ) 


then the system (11) has a periodic solution. Indeed, let Mj be the following compact 
manifold with a boundary 


Since 


d'^V 

dx^ 


x=<5+ln2 


Mi = {xi G M: ( 2 i — 1)(5 + a) ^ Xi ^ 2i{6 + a)}. 
0 and ^ < 0 for X > h + ln 2 , then 


dx 


X=X2 


dV 


dx 


X=Xl 


if (5 + ln 2 ^ xi ^ X 2 - 


Finally, taking into account (12), we obtain that Xj is directed outward to Mj at the boundary 
dMi provided x, = 0 and can be applied. 


3. Conclusion 

Taking into account that the presented approach for studying forced oscillations is rela¬ 
tively simple, and, at the same time, it allows one to use a wide range of ‘building blocks’ 
and types of forces to construct a physical system that has a periodic solution, we hope that 
our result will be useful in applications. 
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